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Summary 


Measurements of constant sagitta, mean gap length versus residual range and photometric 
density versus residual range, have been carried out on 24 stopping protons and 10 z-mesons in 
G5 emulsion. The Fay-Gottstein-Hain p = 0.5 wu (protons) and z = 0.5 mu (a-mesons) scattering 
scheme is calibrated and the conversion factors deduced when calculating the second differences 
from the third and fourth differences respectively. The photometric density is normalised to 
minimum tracks and the absolute normalised photometric coefficient and gap coefficient are 
compared. A linear approximation to the points of photometric density versus residual range in 
a log-log plot is shown to give reasonable values for sufficiently large ranges. An apparatus for 
measuring photometrically ionisation of tracks in nuclear emulsions adapted for the Cooke 
Microscope, is briefly described. 


1. Introduction 


The discovery of a variety of new particles and modes of decay together with the 
use of the new technique of stripped emulsions has made it important to identify 
stopping particles in the emulsions. The possibility of stopping particles has increased, 
and one can follow particles from different events through the stack. The energy 
of the stopping particle can be determined very accurately by measuring the range 
only, provided its mass and charge is known. The charge of the particle can be deter- 
mined by 6-ray counting etc. We restrict ourselves to singly charged particles, however. 
The mass of the particle is another quantity that can be measured. The methods 
most used at present to determine the mass of the stopping particle are based on 
measuring multiple Coulomb scattering, measuring mean gap length versus residual 
range and photometric measurement of ionisation in the track as a function of the 
residual range. Electron sensitive emulsions which begin to saturate at 6 = 0.4 for 
protons, are mostly employed by workers in the field, and thus older methods such 
as blob counting-residual range cannot be used at the end of the track. The method 
of gap counting—residual range can be used but is unreliable. 

The scattering method has the advantage of being independent of any measure- 
ments of ionisation, on the other hand it is usually connected with large statistical 
errors and is sensitive to distortion. All the other methods are based on some estimate 
of the energy losses of the particle in traversing the emulsion. 
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The present investigation has been undertaken to calibrate the Fay-Gottstein-Hain 
scattering schemes and to compare the photometric method and that of mean gap 
length as a measure of ionisation in the determination of the particle mass. Fowler 
and Perkins [17] have used the method of blob and gap counting with great success 
on fast tracks (proton ranges greater than 1 cm) but this method is not applicable 
to the ends of the tracks where the blob density saturates. It was thought that a 
new investigation of the mean gap length near the end of the tracks would be comple- 
mentary to that of Fowler and Perkins. In such an experiment it is important to 
know the exact mass of the particle independently of the measurements otherwise 
extreme signals may be unfairly excluded. The investigation is based on measure- 
ments on 24 protons and 10 z-mesons in an Ilford G5 emulsion, 600 thick, minimum 
blob density 23.9 blobs/100 4 (Sardinia flight no. 21, 1953). The z-mesons are uniquely 
identified from complete 2--e decays. The protons are identified by the consistency 
of the three mass determinations. Two of these are strictly independent. In view 
of the experimental standard deviations one would expect only ~ 1 particle out of 
24 to give a mass value of that of a deuteron at one particular range. Measurements 
are, however, made at 13 different ranges in, for example, the determination of the 
mean gap length and the likelihood of this one track giving an extreme signal at all 
these different points is negligible. The proton sample therefore is a pure one. 

The distribution of track lengths of the observed particles can be seen in Fig. 1. 

All measurements have been performed with a Cooke, Troughton & Simms M 
40180 Nuclear Microscope. 


2. Constant Sagitta method 


(a) Theory 
From the scattering theory it is well known that 
D, 180 z (t \3 
a _ (1) 
t 2 ppc \100 


where ¢ is the cell-size in uw, K the scattering constant, ze the charge, p the momentum 
and fc the velocity of the particle. D, is the absolute mean value of the second dif- 
ferences D{ in uw. No smoothing factor has been taken into account. Fowler’s 
petits method is used. The second and higher order differences can be calculated 
rom 


DP = Yire — 2Yin +Y: 
Df = DP, — DP (2) 
D. =<|D?|) 


y; are the successive readings at points p, of the track at intervals ¢ along the x-axis. 

Scott [1] has given the distribution function of D®, and one can show [2] that 
the mean absolute values of the second, third and fourth differences are related in 
the following manner, 


[DP |) =D, 
(| DP |>=B,= V 20, (3) 
<| Di? |>=2D,. 
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D, measures the deviation of the track from a straight line, D, the deviation from a 
circle with an arbitrary radius and D, the deviation from a curve with linearly varying 
curvature. Thus it is possible to eliminate C-shaped and S-shaped distortion in the 
emulsion by calculating the third and fourth differences respectively. 

The measured mean value D,, consists of true scattering D% and noise 6,. If the 
distributions are normal and true scattering and noise independent, the noise can 
be eliminated according to 

Dn =V(Dr)? — (6n). (4) 


This is approximately true even for higher order differences (S. Biswas et al. [2]). 
The following relation holds between the kinetic energy of a particle (charge ze, 


mass JM) and its range R, 
T =az22" M*-" R*. (5) 


For G5 emulsion with 7’ in MeV, M in proton masses, # in uw the constants are 
a = 0.281 n = 0.568 (3, 4). 
Equations (1) and (5) together with the relativistic expression, 


ppe p 
as =» ( (6) . 
T ~1_Vin# y (6) 
yield 
“ EL K : 
= t!:, ff 
Ds 1800 az?" M'" R"-y (7) 


By choosing, eg ht ete 
fice(l/K) ye" MAR (8) 
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one can maintain the condition D, = const. along the track, i.e. the absolute mean 
value of the second differences of the sagitta is kept constant. The advantages of 
the method are that the number of cells is the largest possible, which means the 
smallest possible statistical error, and a constant signal to noise ratio, so the noise 
is easily eliminated. 

The constant sagitta method has been developed in a number of recent papers 
[5, 6, 7, 8]. 

The varying cell lengths have been calculated by H. Fay, K. Gottstein and K. Hain 
[9] for D,=0.5 uw and different particles, protons, t- and z-mesons (M p — 1836, 
M, = 965, M, = 273 m,), taking into account the variation of the scattering constant 
with cell-size and energy and a 4D — cut-off and assuming a constant noise-level. 
D, =0.5 w has been chosen because it gives a signal to noise ratio about 5. Let 
(Mo, Doo = 9.5 4) be the mass and the absolute mean value of the second differences 
respectively, for the reference particle, and (M, D,) the corresponding values for an 
unknown particle. From equation (7) it can be seen that in the first approximation 
the following relation holds, 


1 
M ) 


iM, \D, iW 
One obtains the relative error in the mass determination, 
o(M)__ 1 _ a(D,) o(D,) 

FT ieens FEe ics 2.01 De (10) 


The relative standard deviation in D, is taken to be 0.75/ VN, where W is the 
number of second differences. This includes statistical errors only. Because of the 
amplication factor 2.31 the statistical standard error in the mass is very large, theo- 
retically of the order of 20% for a proton track of range 5000 su. 


(b) Experimental results 


The scattering measurements have been carried out with an eyepiece filar micro- 
meter [9]. It is necessary to know the noise level of the microscope for the small 
cells employed, down to 10 yw. This has been investigated by measurements on flat 
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Fig. 3. The ratio of the measured third differences to the measured second differences as a function 
of the signal to noise ratio. 


tracks. The tracks have been measured using the small cells, and a large cell, and 
the signal has been eliminated according to the theory of Moliére. The second, third 
and fourth differences have been calculated. The result is shown in Fig. 2. 

The noise indicated includes reading-, stage-, and grain-noise. The errors are standard 
errors derived from the experimental distributions. For the first three points 300 
second differences have been counted, for the last, 100. 

Thus the noise increases with increasing cell size even for small cells. A least square 
fit to the last three points of the second differences yields 6, x t°* in agreement with 
published measurements on the Cooke Microscope [10]. 

Concerning the elimination of noise the following constant values have been 
adopted: 

dy = 0.100 uw, 63 =0.178 uw, 6, =0.315 pw (11) 


corresponding to a cell ¢ = 40 uw. 

The properties of the noise distributions do not seem to have been investigated 
previously. The difference between the properties of the noise distributions and the 
signal distributions is illustrated in Fig. 3. 

When obtaining this figure, 5 flat tracks in different emulsions were measured 
with different cell sizes yielding a total of 2100 second differences. The signal to 


noise ratio was deduced assuming a #2 variation of the signal. 

When calibrating the Fay-Gottstein-Hain scheme the stopping proton and z- 
meson tracks were measured with slight smoothing, 2-3 grains, for small cells (< 40 4) 
and ordinary smoothing, 5 grains, for larger cells. Approximately the same number 
of readings were taken in every track, about 75 for protons and 55 for z-mesons. 
All mean values were calculated with 4 D — cut-off without replacement. Cut-off 
was applied to each of the second, third and fourth differences separately. 

The distributions of the absolute mean values of the differences for the proton 


tracks are shown in Fig. 4. 
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The mean values and relative standard deviations are, 


aD.» = 0.527 fo orl 6 
Oe 

(V2 D;> = 0.545 wo =3.3% 

4D, =0.577 fe o.=3.1%. 

The mean value of D, is in fair agreement with calibration experiments already 
reported [11]. The difference from Friedlander’s, Johnston’s and Merlin’s values, 
0.543, 0.539, 0.549 uw, is due to the different cut-off procedure used. 

The conversion factors between the D,’s calculated, respectively, from the second, 


the third and fourth differences are, 
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Fig. 5. The absolute mean value of the measured second differences as a function of residual range 
for the proton group. 


D, (3rd diff.) _ 

D, Gnd diff) ~ 1-034 + 0.027 
D, (4th diff.) _ 

B, (and diff) ~ 109 + 0.031. 


These values have been calculated without correction for distortion, because it 
has been shown that distortion is negligible in the emulsions. 

In the errors indicated the correlation between the second and following differences 
in the same track has been taken into account. The correlation coefficients are 
003 = 0.66, 05, = 0.57 respectively. 

These results compare favourably with a calibration experiment using other cell 
schemes [12]. There is a definite excess when calculating the signal from the third 
and fourth differences instead of the second and it is necessary to apply the empirical 
conversion factors. 

The corresponding figures for the z-meson tracks are 


i 2 aa 
D, = 0.543 p, (V2-D. = 0.530 u, and <4 Dy) =0.546 w. 


Since they are not of sufficient statistical significance they have not been used when 
calculating the conversion factors. They show, however, that it is necessary to cali- 
brate the 7 =0.5 uw scheme too, and that the deviation from the expected 0.5 yw is 
of the same order as for the proton group. 

The condition D, = const. along the track must be fulfilled if the scattering scheme 
is correct. The variation of the absolute mean value of the measured second differences 
along the track is shown in Fig. 5 for the proton group. 

A weighted least square fit to the points gives D,o< R°°° thus confirming, 
within the limits of experimental error, that D, = sone for the P-scheme. This 
line is shown in the diagram. This point has also been carefully investigated by S. 
Kaneka using u-mesons [13]; he, too, found that the Gottingen scheme is correct. 
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3. Mean gap length versus residual range 
(a) The distribution of gap lengths 


The distribution of gap lengths in tracks of ionising particles in nuclear emulsions 
has been investigated by C. O’Ceallaigh et al. [14, 15, 16] and P. H. Fowler and 
D. H. Perkins [17]. It is shown to be exponential for all measured values of the 
ionisation, i.e. the number of gaps with lengths between x and x+ dz is 


f(x)dx=(1/G) exp (-2/@) dz; (1) 


f(x) is the frequency function of the distribution, G is the mean gap length. 

This distribution is a special case of the Pearson type III distribution. The sample 
distribution of such a population is well known (M. G. Kendall: The Advanced Theory 
of Statistics, vol. 1, p. 244). The sample distribution for the mean value of a sample 
of size N is, if gaps 2 <« are excluded, 


1 (N(a—e)\*! -Xe-oN 
= = = daz. 2 
f(x) dx rm G ec Gi x (2) 
The mean value of this distribution is 
G,=Grte (3) 
and the standard deviation dt 
o = G/N. (4) 


The relations (3) and (4) have been experimentally verified by C. O’Ceallaigh [16]. 
The parameter G, or g = 1/G as used by Fowler and Perkins, seems at present to 
be the best measure of ionisation in the track of a particle. 


(b) The experimental determination of the mean gap length 


Equation (3) makes it possible to measure G in a simple manner. If in a cell ¢ 
with constant ionisation AL, = the total length of all gaps >e« and AN, = the 
number of such gaps, then 

G,=AL,/AN, (5) 


is known and G can be determined. 

Starting from (1) another procedure has been suggested [17]. If N =the number 
of gaps >é¢ per unit length and B =all gaps > 0 per unit length or the blob density 
as usually defined, equation (1) yields 


N =B exp (—¢/G). (6) 


Thus G can be determined by only counting NV and B. As pointed out in the in- 
troduction this method is inapplicable at the ends of the tracks. 

In this paper the first method is used. The choice of ¢ is influenced by several factors. 
The smaller ¢ the more information is obtained from the track, but the measurements 
are very tedious and the subjective errors larger, because one has to judge if a large 
number of small gaps are larger or smaller than ¢. The larger ¢ the less information, 
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Fig. 6. The mean gap length as a function of residual range for protons and z-mesons. 


but the speed is increased. ¢ has been chosen 0.292 yu in order to obtain a sufficient 
number of gaps in cells t = 1000 yw for the proton tracks. 

The track is oriented along the y-axis of the microscope stage and moved forwards 
in the y-direction by the stage micrometer. ¢ is defined by two horizontal lines en- 
graved on the glass plate of an eyepiece filar micrometer ([10] p. 54). The gaps are 
measured with the eyepiece micrometer. In order to facilitate the adjustment of 
the edges of the grains to the horizontal lines the micrometer of the y-axis is equipped 
with a worm gear wheel. 

The mean gap length as a function of residual range is shown in Fig. 6. 

The errors are standard errors calculated from equation (4). It is common practice 
to adopt the same slope for the linear regression line of the particle to be measured 
as for the reference line, in this case the proton line. Using this a weighted linear 


least square fit to the points gives, 


log G=By +B, log R (7) 
for protons By = — 1.1534 B, = 0.2587 
for a-mesons By = — 0.9340 6, = 0.2587. 


These two lines are shown in the diagram. 

For points with equal ionisation, i.e. the same G, the ratio of the masses of the 
particles is the same as the ratio of the residual ranges. If (M, &) is the mass and 
residual range for an unknown particle and (M,, Ry) the corresponding values for 
the reference particle, we have, 

M tithe 
My Ky (8) 
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Fig. 7. The photo-electric apparatus. 


. Movable glass plate 


CO OID OP Ww bo 


. Micrometer screw for the slide 


. The brass tube containing photo-multiplier 10. 
. Eyelens of Cooke’s comp. 8 x eyepiece I i 
. Objective, Cooke’s 5 x 12. 
. Slide for the slit 13. 

Spring 14, 


. Field lens of Cooke’s comp. 8 x eyepiece 15. 


. Brass tube screwed into the microscope 16. 


Slit 

Movable holder with mirror 

Spindle for turning the mirror 

Spindle for turning the glass plate 
Photomultiplier 1 P 21 

Brass tube for the resistances of the 
photo-multiplier 

Electric couplings to the multiplier 


This yields M,/M, = 7.0 which may be considered satisfactory. a 

When comparing with the photometric measurements the gap coefficient g = 1/G@ 
in units mm-’ is used. This is normalised to that of minimum tracks yielding g* = 
9/9o- Jo = 291/mm, measured in the same manner as the coefficient g. 


4. Photometric measurement of ionisation 


(a) The photo-electric apparatus 


An apparatus for the photometric measurement of ionisation has been built. It is 
of the same type as that described by 8. von Friesen and K. Kristiansson [18, 19, 20] 


but adapted for a Cooke microscope. The construction is shown in Fig. 7 


. 


The photo-multiplier is directly coupled to a Speedomax G recorder. With an 
anode voltage of 750 V over the multiplier and a current of 150 «4 Amps suitable 


readings are obtained on the recorder. 
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Fig. 8. Correction factor for the depth of the track in the emulsion. 
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For the measurements a slit 0.3 x 3.0 mm corresponding to 2.9 x 29.0 w in the 
object plane has been used. The turning of the plane parallel glass plate corresponds 
to a displacement of 5 u of the slit on both sides of the track. The paper speed of 
the recorder is 12”/hour. 

If Sj is the scale deflection of the recorder with the slit on one side of the track, 
Si with the slit on the other side and S, with the slit on the track, the ratio 
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Fig. 10. Normalised photometric density as a function of residual range. 
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is measured. The S; are read in 30 yu intervals along the track and the mean value 
S is calculated in successive cells of 600 « for protons and 300 wu for z-mesons. S is 
called the photometric density. 

It is very important to correct for the depth of a track in the emulsion. The correc- 
tion factor is the factor by which a reading should be multiplied in order to be 
normalised to the middle of the emulsion. The correction factor has been deduced 
from the first 600 u of the proton tracks, which are situated at different depths. It 
is shown in Fig. 8. 

It is in fair agreement with the factor deduced by B. Waldeskog [21], but not of 
the same general appearance as that by K. Kristiansson [22]. 

It is highly desirable to normalise the photometric density to that of minimum 
tracks in order to be able to compare the values obtained from different emulsions 
with different degrees of physical development and to compare the values obtained 
by means of other methods of measuring ionisation. This point has been investigated. 

The distribution of 180 photometric densities along two minimum tracks from 
two negatron—positron pairs situated in the middle of the emulsion can be seen in 
Fig. 9. In the same figure distribution of 180 random readings from the same depth 
and in the vicinity of the tracks is shown. 

The mean values are 1.40% and 0.06% respectively, the last value confirming 
that the mean value of random readings equals 0, as it should. The results show 
that it is possible to separate minimum tracks from the background with the slit 
dimensions used. If S, is the photometric density of the minimum tracks, the coef- 


ficient S* = 8/8, is an absolute measure of ionisation in the same sense as g*, normalised 
grain density etc. 
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Fig. 11. The photometric density as a function of residual range in a log-log plot. 


(b) The treatment of the measurements 


The usual way to treat the measurements is to plot the photometric density S 
against residual range R. This is done in Fig. 10. The errors are standard errors derived 
from the spread of the experimental points. 

A parabolic least square curve is fitted to the proton points giving a relation, 


S=fo +p, R + By R* (2) 
with By =12.32 By = —0.661 B,=0.0296 


S in per cent and #& in mm. 

The same calculation for the first five meson points yields 6, = — 3.815 f, = 0.3772. 
The curve is normalised to the same f, as the proton curve. By comparing the residual 
ranges of S = 10.0 in the middle of the measured points one obtains M,/M, = 6.7. 
This may be considered satisfactory and shows that the apparatus is working well. 
Only the first five points of the z-meson curve are used in order to obtain sufficient 
accuracy with a second degree expression. K. Kristiansson has shown [19] that a 
third degree expression is required for L-meson tracks about 1500 yw in length. 

The above treatment is in fact limited to proton tracks less than 10000 uw and L- 
meson tracks less than 2000 u in length, owing to the labour involved in calculating 
higher degree approximations. 

In order to extend the applications of the photometric method and to bring it 
closer to other treatments of ionisation measurements the following procedure is 


suggested. 
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Fig. 12. The relation between the normalised photometric and gap coefficients and residual 
range for protons in G5 emulsion. 


S is plotted versus R in a log-log diagram (no normalisation of the values from 
different particles). Fig. 11 shows that the points are very well approximated by a 
straight line, for 7-mesons with a range > 700 4 and for protons with a range > 5000 uw. 

From a weighted least square fit one obtains a curve 


log S =f, + B, log R (3) 
for protons By = 1.8243 Bp, = — 0.2241 
for 7-mesons Bo = 1.6620 B, = — 0.2353. 


Thus the curves have very nearly the same slope. Adopting the slope of the z-line 
one can determine the M,/M, ratio at 7.4. This is a little too high but it must be 
kept in mind that these calculations are based on very few tracks. The possibility 
of treating the material in this way seems, however, to be established. 


5. Comparison of the normalised ionisation measurements 


The normalised coefficients g* and S* for protons are shown in relation to residual 
range in Fig. 12. The values from the measurements on the z-mesons are used by 
calculating the equivalent proton range assuming a mass of the z-meson of 273 m,. 

For comparison the gap coefficient measured by Fowler and Perkins [17] is shown 
in the diagram. They measured g* by counting blobs and gaps exceeding 4.17 [. 
They obtained identical g* values for different emulsions for ranges longer than 1 em. 
The difference between the g* curves is a real effect, because the g*—as measured 
by the method in this paper—does not differ between two different observers and 
the points are well outside the expected statistical deviation. In order to explain 
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the difference, extensive measurements have been performed in the Bristol B27b 
plate, which has a low degree of development with g, = 186/mm and plateau blob 
density By = 13.5/100 u. Both methods yield the same value of gy. For high ionisation 
g > 1000/mm—there is a definite difference between the measured g values. 

Measurement of the mean gap length gives a g* = 5.3 for 1 cm proton in the Bristol 
plate. Blob and gap counting on the tracks in the Sardinia plate gives g* = 5.7 for 
residual range of 1 em which is in agreement with Fowler and Perkins curve. There 
is thus a shortage of small gaps in the direct measurement of the mean gap length 
compared with blob and gap counting in the region where the blob density begins 
to reach its maximum. That the effect is very pronounced at a range as long as 
1 cm in the Sardinia plate is due to the high degree of development of this plate. 

The normalised photometric coefficient S* is of the same order as the gap coef- 
ficient thus increasing one’s confidence in the normalisation procedure used. The 
curve varies as R-°*. This is somewhat less than the R~°”® variation of the g* 
curve. Fowler and Perkins curve varies as R °**. This indicates that the photometric 
coefficient is less sensitive to the change in ionisation for long ranges than the gap 
coefficient. 


6. Conclusions 


(1) The measured absolute mean values of the second differences for the Fay- 
Gottstein-Hain p= 0.5 uw and a =0.5 uw scattering schemes are p = 0.527 + 0.016 wu 
and a = 0.543 + 0.026 uw; cut-off without replacement being used. (2) The empirical 
conversion factors to be used when calculating the second differences from the third 
and fourth differences are 1.034 + 0.027 and 1.095 + 0.031 respectively. (3) The gap 
coefficient obtained by direct measurement of the mean gap length is smaller than 
the gap coefficient measured by blob and gap counting for g > 1000/mm. (4) The 
relation between the photometric density and residual range is linear in a log-log 
plot for proton ranges exceeding 5000 w. (5) The normalised photometric coefficient 
is less sensitive to the change in ionisation than the gap coefficient for proton ranges 
exceeding 5000 wu. 
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